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We present a systematic study of the variations of the 7 Li NMR properties versus magnetic 
defect concentration ndefcct within the spinel structure of polycrystalline powder samples (nd c f cc t = 
0.21, 0.49, and 0.83 mol%) and a collection of small single crystals (ridefcct = 0.38 mol%) of L1V2O4 
in the temperature range from 0.5 to 4.2 K. We also report static magnetization measurements 
and ac magnetic susceptibility measurements at 14 MHz on the samples at low temperatures. Both 
the 7 Li NMR spectrum and nuclear spin-lattice relaxation rate are inhomogeneous in the presence 
of the magnetic defects. The 7 Li NMR data for the powders are well explained by assuming that 
(i) there is a random distribution of magnetic point defects, (ii) the same heavy Fermi liquid is 
present in the samples containing the magnetic defects as in magnetically pure LiV2 04, and (iii) the 
influences of the magnetic defects and of the Fermi liquid on the magnetization and NMR properties 
are separable. In the single crystals, somewhat different behaviors are observed, possibly due to a 
modification of the heavy Fermi liquid, to a lack of separability of the relaxation effects due to the 
Fermi liquid and the magnetic defects, to non-Fermi liquid behavior of the conduction electrons, 
and/or to quantum fluctuations of finite-size magnetic defects (magnetic droplets). Remarkably, the 
magnetic defects in the powder samples show evidence of spin freezing below T ~ 1.0 K, whereas 
in the single crystals with similar magnetic defect concentration no spin freezing was found down 
to T = 0.5 K. Thus different types of magnetic defects and/or interactions between them appear 
to arise in the powders versus the crystals, possibly due to the substantially different synthesis 
conditions of the powders and crystals. 



I. INTRODUCTION 

LiV204 is a rare ci-electron heavy fermion system at 
low temperatures T < 10 K. 1 The low temperature lin- 
ear electronic specific heat coefficient 7 (0.42 J/mol K 2 ) 
and Pauli magnetic susceptibility \o (~ 0.01 cm 3 /mol) 
are 180 and 310 times those of a free electron gas, respec- 
tively, assuming each vanadium atom contributes 1.5 free 
electrons. The Wilson ratio Rw, which is the ratio of the 
enhancement factors of xo and 7, is equal to 1.7, typi- 
cal for a heavy fermion system^ Heavy fermion behavior 
was further confirmed by electrical resistivity measure- 
ments which show a T 2 dependence below 2 K with a 
large coefficient A = 2.2 [ifl cm/K 2 ^ The A and 7 
values approximately follow the Kadowaki- Woods rela- 
tion, A/7 2 = 1.0 x 10~ 5 fi cm(mol K/J) 2 , which holds 
for a variety of heavy fermion systems £ Despite continu- 
ous theoretical work, a detailed explanation of the heavy 
fermion behaviors in LiV2 04 remains a challenge £££ 

7 Li nuclear magnetic resonance (NMR) was an im- 
portant local probe in establishing the low temperature 
heavy fermion behavior in magnetically pure samples of 
LiV204ii^ The low temperature 7 Li nuclear spin-lattice 
relaxation rate 1/Ti follows a Korringa relation X/T\ on 
T, with a coefficient l/T x T = 2.2 s^K" 1 , which is 6000 
times larger than in the non-heavy fermion isostructural 
superconducting 10 - compound IAT'^Oa^ The Korringa 
ratio re = [ATrk B -/ 2 / (h-y 2 )]K 2 TiT, where K is the Knight 
shift, 7 n and 7 e are the gyromagnetic ratios of the 7 Li 
nuclear spin and the conduction electron spin, respec- 
tively, is equal to 0.7, which is close to the value of unity 
expected for a free electron gas. 

Recently, we found that the low temperature 7 Li 



NMR properties of polycrystalline L1V2O4 are very sen- 
sitive to the presence of a small concentration of mag- 
netic defects (ridefcct = 0.73 mol%) within the spinel 
struct ure i 12 i 13 In a sample containing a negligible con- 
centration of magnetic defects, the longitudinal com- 
ponent of the global bulk 7 Li nuclear magnetization 
M (t) after time delay t following a sequence of satu- 
ration pulses showed a single exponential recovery 1 — 
M(t)/M (00) = exp[-(t/Ti)], where l/T x was propor- 
tional to T as noted above. However, in the sample with 
^defect = 0.73 mol%, the M(i) showed a stretched expo- 
nential recovery 1 - M(t)/M (00) = exp[-(t/T*)P], with 
the characteristic relaxation rate 1/T* showing a peak 
at T w 0.7 K. Here (3 is the stretching exponent with, 
in general, < j3 < 1. There was also a clear differ- 
ence in the 7 Li NMR spectrum in these two samples. At 
low temperatures T < 4.2 K, the magnetically pure sam- 
ple had a narrow spectrum with an almost temperature 
independent width (full width at half maximum peak in- 
tensity FWHM ~ 20 kHz). In contrast, a strong tem- 
perature dependent inhomogeneous broadening (FWHM 
~ 100 kHz at T < 4.2 K) was observed in the sample 
with ridefect = 0.73 mol%. 

In order to further clarify the nature of the magnetic 
defects and their effect on the heavy fermion properties 
of L1V2O4, we report herein 7 Li NMR studies on L1V2O4 
versus magnetic defect concentration. Three polycrys- 
talline samples and a collection of single crystals are stud- 
ied. The powder samples are labeled as 6b, 7a, and 6a, 
with refect = 0.21 mol%, 0.49 mol%, and 0.83 mol%, re- 
spectively. The single crystal sample is labeled as sample 
1 with ^defect = 0.38 mol%. We determined the magnetic 
defect concentrations from static magnetization measure- 
ments in the temperature range 1.8-5 K and applied mag- 
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netic field range 0-5.5 T^z Furthermore, to study the 
possible spin freezing of the magnetic defects at low tem- 
peratures, we measured the ac magnetic susceptibility at 
14 MHz from 0.5 to 6 K of the single crystals and of the 
powder sample 6a with rtdcfcct = 0.83 mol% using the 
tunnel-diode resonator technique?^ 

The temperature dependences of the 7 Li nuclear spin- 
lattice relaxation rates in our polycrystalline samples are 
similar to that of sample 3-3-a2 (ridefcct = 0.73 mol%) 
that we studied in Ref. which showed a peak in 
1/T*(T) at about 1 K. However, we find a qualita- 
tive difference in the temperature dependence of l/Tj* in 
the collection of single crystals, which instead decreases 
monotonically with decreasing temperature from 4.2 K 
down to 0.5 K. We include two important aspects into 
the analysis of the NMR data. First, we consider the 
effect of a distribution of local fields due to different po- 
sitions of the 7 Li nuclei relative to their nearby magnetic 
defects. For the polycrystalline samples, this approach is 
quantitatively consistent with the inhomogeneous broad- 
ening of the spectrum and the nonexponential relaxation 
behavior. In the single crystals, this purely geometric ori- 
gin for the nuclear relaxation fails to explain the observed 
behavior at T < 1.3 K. We then extend our analysis to 
take into account a possible size distribution of postu- 
lated magnetic defects of finite size (magnetic droplets). 
We speculate that the differences between the natures 
and interactions of the magnetic defects/droplets in the 
powders versus the crystals arise from the very different 
synthesis conditions and procedures of the powders and 
crystals. 

The paper is organized as follows. Experimental de- 
tails are given in Sec. [TTJ In Sec. IIII1 we report the ex- 
perimental results of the magnetization, ac susceptibility, 
7 Li NMR spectra, and 7 Li nuclear spin-lattice relaxation 
rate measurements. In Sec. IIV1 we analyze the NMR re- 
sults. In Sec. El we summarize the main conclusions of 
the paper. 



highly sensitive self-resonating LC circuit where losses 
are compensated by a tunnel diode that has a region 
of negative differential resistance in its I-V character- 
istic. The resonant frequency of an empty coil fo = 
1/(27tV LC) changes when a sample is placed in the coil. 
The shift of the resonant frequency, A/ = / (T, H ) — fo 
is directly related to the dimensionless volume ac suscep- 
tibility Xac (T, H) of the sample via 15 



where V s is the sample volume and V c is the coil vol- 
ume. The volume magnetization is the magnetic mo- 
ment per unit volume of the sample, with Gaussian units 
G cm 3 /cm 3 = G. The volume susceptibility is the vol- 
ume magnetization divided by field, which is then di- 
mensionless. The optimized and thermally stabilized cir- 
cuit resonates at 14 MHz with a stability of 0.05 Hz over 
hours!** The resonator was mounted in a 3 He cryostat 
with a temperature range 0.5-150 K. A static external 
field up to 90 kOe can be applied to study field-dependent 
properties. 

7 Li NMR measurements were performed utilizing a 
phase-coherent pulse spectrometer at applied magnetic 
fields H = 1.06, 1.68 and 3.0 T and in the tempera- 
ture range 0.5-4.2 K. Measurements above 1.5 K were 
performed with a 4 He bath cryostat and measurements 
below 1.5 K with a Janis 3 He cryostat. The typical tt/2 
pulse length was 3 jits. The 7 Li NMR spectra for nar- 
row lines (FWHM < 100 kHz) were measured by Fourier 
transform of half the Hahn echo signals, while for wider 
lines, the spectra were measured by integrating the echo 
area as a function of the applied magnetic field at a fixed 
frequency of rf pulses. Nuclear spin-lattice relaxation 
rates were measured by monitoring the recovery of the 
spin echo height using the standard saturation-recovery 
pulse sequence. 



II. EXPERIMENTAL DETAILS 

Polycrystalline LiV2 04 samples were prepared using 
conventional solid state reaction at temperatures up to 
700 °C. The starting materials were V 2 3 (99.99%, MV 
Labs), V 2 5 (99.99%, MV Labs), and Li 2 C0 3 (99.999%, 
Alfa Aeser). Details of the sample synthesis procedure 
can be found in Ref. [U. The typical size of the poly- 
crystalline grains is in the range of 1-10 /im^ as deter- 
mined from scanning electron microscope (SEM) micro- 
graphs. Single crystals were grown at 950-1040 °C using 
a self-flux technique^ 7 - The flux consisted of a mixture 
of Li3V04 and LiV 2 04. The typical size of the crystals 
used in the present work is 0.2 mm. Static magneti- 
zation measurements were performed using a Quantum 
Design SQUID magnetometer in the temperature range 
1.8-350 K and applied magnetic field range 0-5.5 T. 

The ac magnetic susceptibility was measured using a 



III. RESULTS 
A. Magnetic Defect Concentrations 

The magnetic defect concentrations of the samples 
were determined from the low temperature (1.8 K < T < 
5 K) magnetization M versus applied magnetic field H 
isotherms. 14 Figures QJa), []Jb), [ljc), andQJd) show the 
M(H) isotherms at different temperatures for samples 
6b, 1, 7a, and 6a, respectively. The magnetic defect con- 
centration ^defect and spin value S of the magnetic defects 
in each sample are determined by fitting the equation 

M(H,T) = xoH + n dcicct N A gfi B SBs(x) 

= xoH + M do[cct (H,T) (2) 

to all the M(H, T) isotherm data for each sample at 
T < 5 Kfii In Eq. (O, n-dcfcct is the magnetic defect 
concentration in dimensionless mole fraction units, \o 
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FIG. 1: (a)-(d) Magnetization M versus applied magnetic 
field H isotherms at different temperatures T for powder and 
crystals samples of L1V2O4. (e) The magnetic defect contri- 
butions A/defect = M — XoH to the data in panels (a)-(d) 
versus H/(T — 9). The xo and 9 values are listed in Table U 
The solid lines are plots of the second term in Eq. ((2j), Mdefect, 
versus H/(T — 9) with values of ridefect, 9, and S given in Ta- 
ble m 



is a field- and temperature-independent contribution to 
the molar susceptibility at low temperatures T < 5 K, 
Na is Avogadro's number, g is the (powder-averaged) 
spectroscopic splitting factor (g-factor) for the defect 
spins, Bs{x) is the Brillouin function for spin S, and 
x = g[iftSH/[k&{T — 9)] where fee is Boltzmann's con- 
stant. We have replaced T in the usual Brillouin func- 
tion by (T — 9) in order to take into account weak 
interactions between the magnetic defects. \o, S, 9, 
and ridefect are free parameters in the fit, whereas the 
g value is fixed to be equal to 2 during the fitpii The 
best fit parameters are listed in Table Q] Figure [lje) 
shows the magnetic defect magnetization contributions 
M dcicct (H,T) = M(H,T) - xoH versus H/(T - 9) for 
the four samples. All the data points in Figs. [lja), [ljb), 
[He), and [ljd) fall onto a universal curve in Fig. [He) 
for each sample, respectively, as described by the second 
term in Eq. @, thus confirming the consistency of the 
fits. The fitted functions M de{ect (H,T) in Eq. © for the 
four samples are plotted versus H/(T — 9) as the solid 
curves in Fig. [lje) and show excellent agreement with 
the data. 



TABLE I: Best fit values of the magnetic defect concentra- 
tion ridefect, the spin value S, the intrinsic susceptibility xo, 
and the effective Weiss temperature 9 for powder samples 
6b, 7a, and 6a and crystal sample 1 of LiV204 obtained by 
fitting Eq. (0 to the low temperature (1.8 K < T < 5 K) 
magnetization versus field isotherms in Figs. [Ha)-[TJd) with 
< H < 5.5 T. 



Sample 


"-defect (mol%) 


S 


Xo (cm 3 /mol) 


9(K) 


6b 


0.21(1) 


3.6(2) 


0.0104(1) 


-0.75(14) 


7a 


0.49(1) 


3.5(1) 


0.0108(1) 


-0.57(6) 


6a 


0.83(3) 


3.9(1) 


0.0122(2) 


-0.64(10) 


1 


0.38(1) 


3.3(1) 


0.01186(4) 


-0.43(6) 



Several features of the data in Tableware important. 
The low-temperature field-independent (up to 5.5 T) sus- 
ceptibilities Xo °f a U f° ur samples are the same to within 
about 10%, even though the magnetic defect concentra- 
tions change by a factor of 4, and are about the same 
as in magnetically pure LiX^O^^ This agreement sug- 
gests that the heavy Fermi liquid in magnetically pure 
LiV2C>4 survives in the presence of the magnetic defects. 
Second, the spins S of the magnetic defects should be 
considered as average values, and these values are large, 
ranging from 3.3 to 3.9. That the magnetic defects have 
large spins is obvious from the data in Fig. [IJe) because 
the magnetic defect magnetizations are nearly saturated 
at relatively low fields of only ~ 2 T; spins 1/2 would 
not saturate even at our maximum field of 5.5 T. It is 
difficult to understand how such large spin values could 
arise from point defects in the crystal structure. In that 
case one might expect the magnetic defect spins to be 
much smaller and similar to those of V +4 (S = 1/2) or 
V +3 (S — 1). The large spins of the magnetic defects thus 
suggest that these spins may be associated with extended 
objects that we call "magnetic droplets" in Sec. lIVCl be- 
low instead of being associated with point-like local mag- 
netic moments as in the usual picture. Third, the Weiss 
temperatures 9 for all the samples are rather small, and 
indicate that the average interaction energy between the 
magnetic defects is also small and of order 1 K. Finally, 
from low-temperature magnetization measurements on 
many polycrystalline and single crystal samples that we 
have carried out in addition to those described here, the 
magnetic defect concentrations found do not exceed the 
largest value listed in Table Q] of 0.83 mol%. 



B. ac Magnetic Susceptibility at 14 MHz 

The ac magnetic susceptibility, Xac = dM/dH, is an 
important parameter directly related to the electronic 
spin dynamics. It is very sensitive to collective behavior 
such as spin freezing and a transition to the glassy state. 
Figure [5] shows Axac = Xac (T) - Xac (4.8 K) versus tem- 
perature T at various values of the external magnetic 
field for powder sample 6a with n.d c f cct = 0.83 mol%. 
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FIG. 2: (color online.) A^ac, the change of ac magnetic 
susceptibility Xac at 14 MHz relative to its value at 4.8 K, 
versus temperature T for L1V2O4 powder sample 6a with 
rtdcfcct = 0.83 mol% at several values of the magnetic field 
(indicated in the legend). Xac decreases with increasing mag- 
netic field. 



Each curve corresponds to a magnetic field listed in the 
legend and the curves from top to bottom correspond to 
increasing magnetic field. We note that the change of 
the magnetic moment amplitude of the measured sam- 
ple on decreasing the temperature from 1.1 K to 0.5 K in 
zero field corresponds to a change in magnetic moment of 
only about 5 x 10 -10 G cm 3 , which cannot be resolved by 
a conventional SQUID magnetometer for the same size 
(~ 0.3 mm 3 ) sample. 

At zero static applied field, there is an obvious peak 
in Axac at about 1.1 K in Fig. [2] that is most likely in- 
dicative of a collective freezing of the magnetic moments. 
The field dependence of the magnetic susceptibility is 
characteristic of a spin glass system where spin random- 
ness is suppressed by the uniaxial field and the peak in 
Xac associated with spin freezing is suppressed because 
the magnetic moments are closer to saturation. This re- 
sult suggests collective freezing behavior of the magnetic 
defects in the L1V2O4 powder sample in zero field. 

For our sample 1 consisting of a collection of single 
crystals with overall ndcfcct = 0.38 mol%, the situation is 
quite different. We cannot measure the spin susceptibil- 
ity because the diamagnetic orbital susceptibility Xac.skin 
arising from skin depth effects dominates it. The skin 
depth 8 can be calculated fromi£ 



6 = 



504 



meters. 



(3) 



where K m is relative permeability, er is the conductivity 
in fi m , and v is the applied frequency in Hz. Set- 
ting ifm = 1, cr = 5 X 10 6 Q- 1 m- 1 (a value at 1.8 K 
in Ref. fl7t) . and v = 14 MHz, we obtain S w 0.06 mm, 
significantly smaller than the size of a crystal in sam- 
ple 1. Thus we expect that the Xac,skin contribution to 



3 



1— 

CO 



1 1 1 1 1 1 1 1 1 1 


1 1 1 


single crystal 1 . 




H = / 


- 


— ^1 


, 1 , 



1 2 3 4 5 6 7 
T(K) 

FIG. 3: Axac, the change of ac magnetic susceptibility Xac 
relative to its value at 0.5 K, versus temperature T in single 
crystal L1V2O4 (sample 1) in zero applied magnetic field and 
at a frequency of 14 MHz. 



Xac is significant and its effect increases with decreas- 
ing temperature as the resistivity decreases monotoni- 
cally with decreasing temperature; 3 "! 4 - Figure [3] shows the 
Axac = Xac(T) — Xac(0.5 K) versus temperature T from 
0.5 to 6 K. Since the static susceptibility of various sam- 
ples is nearly T-independent or increases with decreasing 
T over this T range, the decrease in Axac with decreasing 
T in Fig. [3] indicates that Xac, skin (T) dominates the Xac 
response there. Furthermore, we see no evidence for a col- 
lective spin freezing for this sample, and we did not find 
any field dependence up to an applied field of 10 kOe (not 
shown). Thus, the measurement of the ac susceptibility 
at 14 MHz for our single crystals does not yield useful 
information for understanding the magnetic response of 
the magnetic defects in these crystals. Xac(H,T) mea- 
surements at much lower frequencies are called for. 



C. 7 Li NMR Line Width 

The 7 Li NMR absorption line width is related to the lo- 
cal static magnetic field distribution. It becomes broader 
with increasing concentrations of magnetic defects. Fig- 
ure [4] shows the absorption lines of the four samples at 
temperature T = 4.2 K and H = 1.06 T. Although the 
7 Li nuclei have spin / = 3/2, both first and second or- 
der nuclear quadrupole broadening due to a structural 
distortion can be ruled out since we observe no satel- 
lite peaks or shortening of 7r/2 pulse length as compared 
to the magnetically pure L1V2O4 samplei 19 ' 20 The line 
width is significantly larger than the intrinsic width for 
an individual 7 Li nuclear spin, indicating an inhomoge- 
neous magnetic broadening of the line. The intrinsic line 
width is of the order of I/T2 ~ 5 kHz, where T2 is the 
nuclear spin-spin relaxation time and is almost indepen- 
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FIG. 4: The 7 Li NMR absorption versus rf frequency v at 
temperature T = 4.2 K and applied magnetic field H = 1.06 T 
in the four LiV2C>4 samples. The frequency vo = 17.6 MHz. 



dent of the defect concentration and temperature below 
4.2 K. Figure [5] displays the temperature dependences of 
the full width at half maximum peak intensity (FWHM) 
of the spectra for the four samples. 

The broadening of the 7 Li NMR line has three con- 
tributions. The first contribution comes from the nu- 
clear 7 Li- 51 V and 7 Li- 7 Li dipolar interactions. This con- 
tribution can be estimated using the Van Vleck sec- 
ond moment (Auj 2 )^- A second broadening comes from 
the macroscopic field inhomogeneity due to a distribu- 
tion of the demagnetization factors and a distribution 
of magnetic fields due to neighboring powder grains. 
This contribution is proportional to the magnetization 
of the sample and the resulting root mean square devi- 
ation of 7 Li NMR resonance frequencies can be written 
as BMpNju/Zn, where M is the molar susceptibility, 
Pn is the density of L1V2O4 formula units in the sam- 
ple, 7Li the gyromagnetic ratio of 7 Li nuclei, and B a 
dimensionless factor. B is estimated to be 1.43 for a 
close packed powder sample with ellipsoidal shapes^ A 
third broadening contribution comes from inhomogene- 
ity due to the presence of magnetic defects within the 
sample. An estimate for this contribution is not possi- 
ble without a model of the nature of the defects and the 
types of interactions between the defects and nearby 7 Li 
nuclear spins. However, the presence of this contribution 
can be inferred by comparing the experimental FWHM 
values and the values expected when including only the 
first two contributions, as follows. 

The FWHM resulting from the first two contributions 
can be calculated within a Gaussian approximation by 



FWHM a = 2.35v / (Aw 2 )/(2tt) 2 + (BM p Nllji l2Ttf (4) 

with B = 1.43, and (Aw 2 ) 1 / 2 /2 7 r = 2.7 kHz. 23 M is cal- 
culated from Eq. ([2]) using the parameter values listed in 
Table H The FWHM a calculated from Eq. © is plotted 
as the dashed lines in Fig. \5\ It is clear that Eq. can- 
not account for the observed broadening of the lines, so 
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FIG. 5: Temperature T dependence of full width at half max- 
imum peak intensity (FWHM) of the 7 Li NMR spectrum un- 
der external magnetic field H = 1.06 T in the four LiV2 04 
samples. The symbols are experimental results. The dashed 
lines are plots of Eq. Q (with B = 1.43) that takes into ac- 
count the contributions due to powder broadening and nuclear 
dipole-dipole interactions, but does not take into account lo- 
cal field inhomogeneity due to the magnetic defects. The solid 
lines are fits by Eq. (0, which also takes into account the lo- 
cal field inhomogeneity. The fitted solid lines from bottom to 
top are for samples with ridefoct = 0.21, 0.38 (crystals), 0.49, 
and 0.83 mol%, respectively. 



a local magnetic field inhomogeneity due to the presence 
of the magnetic defects must be present in the samples. 
We will return to this issue in Sec. IIV B 11 



D. Nuclear Spin-Lattice Relaxation Rates 

The longitudinal 7 Li nuclear spin relaxation versus 
time M(t) exhibits an increasingly nonexponential be- 
havior with increasing concentration of magnetic defects 
or decreasing temperature. Figure [5] shows the recoveries 
of M(t) following a saturation sequence for the four sam- 
ples at different temperatures. The recovery data can be 
described by a stretched exponential function 



M(t) 
M(oo) 



exp[-(t/T*; 



(5) 



The solid curves in Fig. [75] are best fits to the data by 
Eq. ©. The best fit values of l/Tf(T) and (3(T) are 
shown in Figs. [7] and [8] for powder and single crystal 
samples, respectively. 

The temperature dependence of 1/T-j* is quite dif- 
ferent in the powder and single crystal samples. A 
peak is observed in 1/T*(T) for the powder samples 
6a (ndefect = 0.83 mol%, T pcak ~ 1.0 K) and 7a (n dcfcct = 
0.49 mol%, T pcak S3 0.6-0.7 K). In the powder sam- 
ple 6b with the smallest magnetic defect concentration 
("defect = 0.21 mol%), 1/T* starts to increase at the low- 
est experimental temperatures and might exhibit a peak 
with further decreasing temperature. The peak positions 
in sample 6a for H — 1.06 and 1.68 T are almost the 
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FIG. 6: Recovery of 7 Li nuclear magnetization M(t) after 
time delay t following a sequence of saturation pulses. Note 
that M here is different from the electronic spin magneti- 
zation in Fig. [TJ The data points were obtained in applied 
magnetic field H = 1.06 T at the indicated temperatures and 
with rf frequency v = 17.6 MHz for L1V2O4 samples with (a) 
0.21 mol%, (b) 0.38 mol% (single crystals), (c) 0.49 mol%, 
and (d) 0.83 mol% magnetic defects. The solid curves are fits 
to the data by Eq. ([5]). 



same as the peak position in XacC? 1 ) for this sample at 
H = in Fig. H We conclude that the peaks in 1/T-* 
originate from the spin freezing of the magnetic defects. 
In the crystal sample, 1/T*(T) in Fig. [S] decreases mono- 
tonically with decreasing temperature with a 1/T* value 
at 0.5 K much smaller than in the powder samples, and 
there is no sign of spin freezing. 

Before ending this subsection, we comment about the 
effect of inhomogeneous broadening on the relaxation 
measurements. Because of the increasing inhomogeneous 
broadening with decreasing temperature, some of the 
7 Li nuclei may be shifted out of the NMR spectrom- 
eter response window (A/ ~ 200 kHz) and excluded 
from the relaxation measurements. The number of ob- 
served 7 Li nuclei can be estimated from the product of 
fully recovered echo height M(oo) and the temperature, 
which is proportional to the nuclear Curie constant C 
in the Curie law for M(oo) = C/T. These data are 
shown for H = 1.06 T versus temperature T in Fig. [5] 
For powder samples 6b (rtdefoct = 0.21 mol%) and 7a 
(^defect = 0.49 mol%), the decrease of M(oo)T is less 
than 10% when the temperature decreases from 4.2 K 
to the lowest temperature (« 0.5 K). In contrast, for 
sample 6a (ndefect = 0.83 mol%), M(oo)T starts to de- 
crease below T ps 3.5 K and at the lowest temperature 
(T fa 0.5 K), M(oo)T is about 50% of that at 4.2 K. As 
we will show below, the nuclei at the wings of the spec- 
trum have an average relaxation rate larger than those 
at the center of the spectrum. Exclusion of those nu- 
clei in sample 6a can thus result in a smaller measured 
relaxation rate in that sample. 

In the single crystals, the normalized signal inten- 
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FIG. 7: 1/T* and f) versus temperature T of L1V2O4 obtained 
by fitting data as in Figs. EJa), (c), and (d) by Eq. ©, of (a) 
powder sample 6b with ridefect = 0.21 mol% at external mag- 
netic fields H — 1.06, 1.68, and 3.0 T, (b) powder sample 7a 
with ndefect = 0.49 mol% at H = 1.06, 1.68 T, and (c) powder 
sample 6a with nd c f cc t = 0.83 mol% at H = 1.06 and 1.68 T. 



sity M(oo)T also decreases with decreasing temperature. 
Since the line width in the crystals is less than in powder 
sample 7a (see Fig.[S]), where no significant signal loss is 
observed, we attribute the signal loss to the effect of rf 
field skin depth. Here, only the 7 Li nuclear spins within 
the skin depth contribute to the NMR signal. Setting 
K m = 1, or = 5 x 10 6 ft^m- 1 (the value of a at 1.8 K in 
Ref.Qj]), and v = 17.6 MHz, Eq. © gives S = 0.054 mm, 
which is less than the typical size (0.2 mm) of the crys- 
tals. However, there is an unexplained kink in the data 
for the crystals at T ps 1.4 K in both Figs. [5] and [U 
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FIG. 8: 1/T* and f3 versus temperature T of the LiV2 04 crys- 
tal sample 1 with ndcfcct = 0.38 mol% in external magnetic 
fields H = 1.06 and 1.68 T, obtained by fitting data as in 
Fig. Ob) by Eq. ©. 
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FIG. 9: The fully recovered echo intensity M(oo), which is 
the total equilibrium nuclear magnetization, times temper- 
ature T versus T in the nuclear spin-lattice relaxation rate 
measurements of the four LiV2 04 samples in an applied field 
H = 1.06 T. 
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FIG. 10: Recovery at 4.2 K of a "hole" in the absorp- 
tion spectrum of LiV2 04 produced by a weak 7r/2 pulse 
with pulse length of 56 its at delay = in sample (a) 6a 
(n dcfcct = 0.83 mol%) and (b) 6b (n dc f cc t = 0.21 mol%). The 
applied magnetic field H is 1.06 T and the center frequency is 
17.6 MHz. The delay times after which the spectra were mea- 
sured by two strong rf pulses are given in the figures. Note 
the different abscissa scales in (a) and (b). 



E. Relaxation at Different Positions in the Spectra 

The observation of a stretched exponential relaxation 
behavior indicates the presence of a distribution of nu- 
clear spin-lattice relaxation rates 1 /T\ . In order to study 
the origin of the 1 /7\ distribution, we performed the fol- 
lowing "hole burning" experiment. This experiment ex- 
tends our previous hole burning experiment briefly de- 
scribed in Ref. We also studied the relaxation behav- 
ior at different positions of the NMR absorption line. 

Figures [TOT a) and (b) display the recovery of a "hole" 
in the echo spectrum in applied magnetic field H = 
1.06 T, obtained from Fourier transform of half the Hahn 
echo signal generated by two strong rf pulses following 
a weak ir/2 pulse in samples 6a (ndcfcct = 0.83 mol%) 
and 6b (n^efect = 0.21 mol%), respectively. The weak 
7r/2 pulse has a width of 56 fj,s and most of its power is 
distributed within a narrow frequency window of width 



w 40 kHz. Such a weak 7r/2 pulse only saturates the 
central part of the spectrum. It is clear that the hole 
recovery process does not affect the rest of the line and 
thus spectral diffusion does not occur in our time scale. 
That is, nuclei with different Larmor frequencies are not 
coupled to each other over the NMR relaxation time scale 
of Ti - 100 ms. 

Lack of spectral diffusion as observed above allows us 
to investigate the nuclear spin-lattice relaxation at dif- 
ferent positions of the spectrum. Due to the strong Li 
NMR signal at low temperatures, we were able to study 
the relaxation of 7 Li far out on the wings of the spectrum 
although the signal intensity is much weaker than at the 
peak. Figure [TT] displays the nuclear spin-lattice relax- 
ation curves of powder sample 7a (ndefect = 0.49 mol%) in 
H = 1.68 T with the rf pulse frequency equal to, 400 kHz 
higher than, or 400 kHz lower than, the peak frequency 
of the line. All three recovery curves are nonexponen- 
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FIG. 11: Recovery at 4.2 K of 7 Li longitudinal nuclear mag- 
netization M(t) following a saturation sequence at time t = 
measured at different positions of the spectrum in LIV2O4 
powder sample 7a (ridefect = 0.49 mol%) under external mag- 
netic field H = 1.68 T. The recovery curves, which are non- 
exponential, were measured with rf pulse frequency (□) equal 
to the peak of the spectrum (27.8 MHz) and at T = 1.56 K, 
(A) 400 kHz lower than the peak and at T — 1.53 K, and 
(o) 400 kHz higher than the peak and at T = 1.77 K, 
respectively. 



tial. It is clear from Fig. [TT] that the nuclei close to the 
peak of the line have an average relaxation rate lower 
than those away from the peak. As will be discussed 
below, the behavior in Fig. [TT]is consistent with an inho- 
mogeneous local magnetic field induced by the magnetic 
defects. It is noted that the temperatures at which the 
three relaxation curves were taken are slightly different. 
However, such small temperature differences should be 
negligible compared to the large difference of relaxation 
rates between these three curves. 



IV. ANALYSIS 



A. Introduction 



We will assume at the outset that the heavy Fermi 
liquid present in magnetically pure L1V2O4 is not af- 
fected by the presence of the magnetic defects. The 
measured 7 Li relaxation versus time following satura- 
tion then arises from two mechanisms. The first is a 
single-exponential relaxation that is the same as in mag- 
netically pure L1V2O4 and that comes from the contact 
interaction of 7 Li nuclear spins with the conduction elec- 
tron spins in the heavy Fermi liquid. The second mech- 
anism is the hyperfine interaction of the nuclei with the 
magnetic defects. The first (homogeneous) mechanism 
gives a relaxation rate described by the Korringa law 
with 1/Ti cx Tj2i The main goal of the present mod- 
eling is to then determine the (second) contribution of 
the magnetic defects to the time-dependent nuclear re- 
laxation and to subsequently interpret what that contri- 
bution means. The separation of these two contributions 
to the magnetic properties, at least above 1.8 K, is sup- 
ported by previous magnetization measurement s 1 ' 1 as 
well as by those in Sec. IIII Al above. The magnetization 
as expressed in Eq. ^ contains a contribution xoH al- 
most independent of the magnetic defect concentration 
(see Table fl]) . This contribution is most likely due to the 
same heavy Fermi liquid that is present in magnetically 
pure L1V2O4 at low temperatures. 

We will see that this separation of the magnetic proper- 
ties into a heavy Fermi liquid part and a magnetic defect 
part can consistently explain our 7 Li NMR measurements 
on our powder samples of L1V2O4. However, as we will 
show in Fig. fl~5Ta) below, our NMR longitudinal mag- 
netization recovery data below ~ 1.3 K for our sample 
of single crystals indicate that the Fermi liquid is modi- 
fied by the presence of magnetic defects and/or that our 
model for the magnetic defects is no longer accurate be- 
low that temperature in our single crystals. This differ- 
entiation between the bulk crystal and powder properties 
indicates that there are differences between the natures 
of the magnetic defects and/or their interactions in single 
crystals as compared to powders, which in turn are likely 
associated in some way with the quite different prepara- 
tion conditions of the two types of samples. 



A model for the microscopic nature of the magnetic 
defects has to be assumed in order to analyze the NMR 
results. We will examine two related possibilities. First, 
in Sec. HVBl we treat the magnetic defects as traditional 
identical localized magnetic moments. The distribution 
of 7 Li nuclear spin-lattice relaxation rates and the inho- 
mogeneous broadening of the 7 Li absorption spectrum 
are then entirely caused by the local field inhomogeneity, 
which arises from a distribution of positions of the 7 Li 
nuclei relative to the magnetic defects. In a refined ver- 
sion of this model in Sec. IIVC1 the magnetic defects are 
assumed to actually be "magnetic droplets" that have a 
distribution of sizes. In this refined approach, the dis- 
tribution of the dynamic properties of the droplets also 
needs to be considered. 



B. Geometric Inhomogeneity 

1. 7 Li NMR Line Width 

First we will analyze the 7 Li NMR line width by consid- 
ering the distribution of distances between nuclear spins 
and point-like magnetic defects within the spinel struc- 
ture. Dilute paramagnetic centers give rise to a broaden- 
ing of the NMR spectrum through inhomogeneous dipo- 
lar and RKKY interactions and in the limit of great di- 
lution the line shape approaches a Lorentzian with full 
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width at half maximum intensity (FWHMb) given by 2 ^ 

8np N 



FWHM b = An dl 



feet ' 



9\/3 



4.5An dclcct SB s (x) MHz, 



(6) 



where pn = 1-44 x 10 22 cm -3 is the number density 
of L1V2O4 formula units, A = 1 for purely dipolar in- 
teractions and A > 1 if the RKKY interaction is also 
important, (S z ) is the thermal average value of mag- 
netic defect spin polarization along the direction of the 
applied magnetic field and is equal to SBs(x) with 
x = gp B SH/[k B (T - 9)] [see Eq. ©]. The line shape 
due to the dilute magnetic defects is Lorentzian 2 ^ while 
the line shape due to the two contributions in Eq. (g]) is 
Gaussian i 21 ' 22 In order to obtain the final FWHM value, 
we convolute a Gaussian distribution with FWHM = 1 
with a Lorentzian distribution that has FWMH = x and 
the same mean value as the Gaussian distribution. We 
find that the FWHM of the convoluted distribution can 
be approximated by (1 -I- x 8 / 5 ) 5 /8 to within 10% for all 
values of x. We estimate the total FWHM by combining 
Eqs. (|4]) and ([6]) according to 

FWHM = (FWHM a t+FWHM b 5)l 

5. 

= {FWHM a t + [A.5An dcicct SB s (x) MHz] ^ j 8 . 

(7) 

Using the values of n dofect , S and 6 in Table U and the 
results for FWHM a in Fig. E3 the FWHM data in Fig. [5] 
for all four samples were simultaneously fitted by Eq. ([7]), 
except for the single crystal data below 1.3 K, where 
the nuclear spin-lattice relaxation rates in Fig. [5] indicate 
a possible screening of the magnetic defects. The only 
fitting parameter was A, and the best fit value was A = 
1.4. The best fit to the data is displayed as the set of 
solid curves in Fig. The high quality of the fit shows 
that the local field inhomogeneity at the 7 Li nuclear sites 
arising from the distribution of distances between the 7 Li 
nuclei and the magnetic defects is an essential contributor 
to the 7 Li NMR line width. 



2. 7 Li Nuclear Spin-Lattice Relaxation 

In the present approach we treat the magnetic defects 
as identical localized magnetic moments. The distribu- 
tion of 7 Li nuclear spin-lattice relaxation rates 1/Ti then 
arises from a distribution of fluctuating local magnetic 
fields at the nuclear sites due to a distribution in the 
positions of the nuclei relative to the magnetic defects. 
Since the relative positions of the 7 Li nuclei with respect 
to the magnetic defects are fixed, the shape of the 7 Li 
1/Ti probability distribution due to the defects should 
be temperature independent. This would give rise to a 
temperature independent (3 valu o 12 i 26 in the stretched ex- 
ponential function in Eq. (JSJ) if there were no additional 
contributions to the 7 Li nuclear spin-lattice relaxation. 



The observed temperature dependences of the stretch- 
ing exponent (3 in the insets of Figs. [7] and [5] are explained 
in this model by the additional Korringa contribution to 
1/Ti that is proportional to the temperature. Since the 
nuclear spin-lattice relaxation rate due to itinerant con- 
duction electrons is assumed to be homogeneous across 
the sample since it results from the contact interaction 
between the nuclear and conduction electron spins, the 
nuclear spin recovery due to the conduction electrons 
alone should be a single exponential. As just discussed, 
the recovery due to the magnetic defects alone should 
be a stretched exponential function with a temperature- 
independent (3. The observed temperature dependent (3 
arises in our model from different temperature depen- 
dences of the Korringa and magnetic defect contributions 
to the nuclear spin-lattice relaxation. Different tempera- 
ture dependences result in different weights of these two 
contributions at different temperatures and accordingly 
different j3 values are seen at different temperatures when 
the total recovery is fitted by a stretched exponential 
function Eq. ([5]). Similarly, it is not appropriate to an- 
alyze the 1/T X *(T) data in Figs. [JJ and [5] in terms of a 
sum of contributions from the heavy Fermi liquid and 
from the local magnetic defects, because their respective 
contributions to 1/T^(T) cannot be deconvoluted. 

To determine the magnetic defect contribution to the 
7 Li nuclear spin dynamics, we first extract from the 
observed 7 Li nuclear spin-lattice relaxation versus time 
M{t) data what the contribution of the magnetic defects 
is, and then derive parameters describing the relaxation 
by the magnetic defects. To accomplish the former goal, 
we write 1- M(t)/M (oo) = p(t) exp(-t/T lK ), where p(t) 
is the contribution to the 7 Li nuclear spin relaxation from 
the magnetic defects and exp(— t/Tix) is the Korringa 
contribution from the heavy Fermi liquid, where we as- 
sume a concentration independent Korringa relaxation 
rate l/Ti K = (2.2 s _1 K _1 )T and the coefficient of T is 
taken to be the value in a magnetically pure sampled 
Then one obtains 



P(t) 



1 



M(i) 



M(oo) 



exp 



t 



T 



IK 



(8) 



Thus p(t) is determined by multiplying the experimen- 
tally observed 1 - M(t)/M(oo) by exp(t/Ti K ). 

We find that the magnetic defect contribution p(t) to 
the 7 Li nuclear spin-lattice relaxation usually follows a 
stretched exponential time dependence in our tempera- 
ture range 0.5 < T < 4.2 K with a temperature- and mag- 
netic defect concentration-independent stretching expo- 
nent /?, as anticipated above, where we find that (3 has 
the specific value (3 = 1/2. Thus we obtain 



p(t)=exp[-(t/I7 d ) 1/a ] 



(9) 



where the new parameter T* d takes the place of T* in 
Eq. ©. Figures H2] and UMsl) show plots of the log- 
arithm of pit") versus i 1 / 2 in external magnetic field 
H = 1.06 T and at different temperatures for powder 
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and single crystal samples, respectively. In powder sam- 
ples 7a (ridefect = 0.49 mol%) and 6a (ndefect = 0.83 
mol%), p(t) can be fitted very well by Eq. (|9|) at all 
temperatures, as shown by the linear fits in Figs. 112( b) 
and (c) , respectively. In powder sample 6b with a smaller 
^defect = 0.21 mol% in Fig. fT27a). p(t) follows root ex- 
ponential behavior for all times t at the higher tempera- 
tures but only at short times at the low temperature of 
0.61 K. We infer in Sec. IIV B~4l below that the deviation 
at longer times is due to the effect of spin diffusion. In the 
crystals, p(t) in Fig. fTBT a) follows root exponential decay 
above 1.3 K but at lower temperature p(t) instead shows 
an unphysical increase at later times. This unphysical 
behavior suggests that Eq. (|8|) overestimates the conduc- 
tion electron contribution to the nuclear spin-lattice re- 
laxation at temperatures below 1.3 K, the separability of 
the relaxation due to the magnetic defects from that due 
to the conduction electrons is no longer appropriate in 
that temperature regime, or the conduction electrons no 
longer form a Fermi liquid. An additional possible reason 
for the unphysical behavior is given in Sec. IIV 01 Resolv- 
ing this issue is an important topic for future research. 

We extract 1/Tj^ versus temperature from the slopes 
of the fitted lines of log[p(i)] versus t 1 / 2 in Figs. [12] and 
WSia) according to Eq. The results are displayed in 
Figs.fTBTb) and[T4lfor the single crystal and powder sam- 
ples, respectively. The 1/T* d versus T in powder sam- 
ples 7a and 6a in Figs. fT4Tb) and (c), respectively, show 
an almost field independent peak, similar to the peaks in 
1/T* versus T in Figs.[7[b) and (c). As discussed above, 
the peaks are attributed to spin freezing of the magnetic 
defects. For the single crystals, we only extract 1/T* d 
values above 1.3 K for reasons discussed above. Here 
1/T* d is nearly constant from 4.2 K down to about 2 K, 
but then shows a decrease upon further decrease in T. 
This behavior is very different from that of the powders. 

The above root exponential relaxation behavior has 
been reported previously in systems where the nuclear 
spin-lattice relaxation rate is proportional to 1/r 6 , where 
r is the distance between a nucleus and a nearby para- 
magnetic center, and no nuclear spin diffusion takes 
placei 27 ' 28 Nuclear spin-lattice relaxation due to fluctua- 
tions of both dipolar and RKKY interactions have such 
1/r 6 dependences. In general, one can write the nuclear 
spin-lattice relaxation rate at the nuclear site r due to a 
nearby paramagnetic center at location R; as 



1 



Zi(r,R,) 



f(0) 



1 



/ |r-Rz 



(10) 



where r — R; is the vector connecting the paramagnetic 
center and nuclear spin, 9 is the angle between r — R; 
and the external magnetic field, f(0) is the angular de- 
pendence of 1/Ti(r, R ; ), / is the average of f(0) over all 
directions, and 



Ci — C - 



n 



(11) 

1 + {^vTlY 

is a parameter proportional to the spectral density of spin 
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FIG. 12: Semilog plot of p(t) in Eq. ((8} versus the square root 
of the delay time t 1 ^ 2 for the LiV2C>4 powder samples (a) 6b 
with ridcfoct = 0.21 mol%, (b) 7a with ndefect = 0.49 mol%, 
and (c) 6a with ndefect = 0.83 mol% at applied magnetic field 
H = 1.06 T and different temperatures. The straight lines 
are best fits of the data by Eq. ©, with parameters 1/T* d 
given in Fig. 1141 



fluctuations (with correlation time r;) at the nuclear Lar- 
mor frequency o)„i21 The defect contribution of the recov- 
ery of the nuclear magnetization towards equilibrium is 
then given as 



p (r, t) = exp 



Ti(r,R,) 



(12) 



where the sum is over all defect sites. We first ignore the 
spatial variation of the microscopic relaxation rate, i.e. 



11 



— I — I — I — I — I — I — I — I — I — I- 
H = 1.06 T 
crystal sample 1 
n , = 0.38 mol% 

defect 




i i_ 



1.5 



t 1 ' 2 (ms 1 ' 2 ) 




^ 0.5 - 



T(K) 

FIG. 13: (a) Semilog plot of the relaxation function pit) in 
Eq. ([8]) versus the square root of the delay time t 1 ^ 2 for the 
LiV204 crystal sample 1 with ndcfcct = 0.38 mol% at applied 
magnetic field H — 1.06 T and at two different temperatures. 
The upturn in the 0.79 K data at large times is unphysical 
(see text). The straight line is a best fit of the 3.92 K data 
by Eq. JS). (b) l/T* d in Eq. © versus temperature T in 
H = 1.06 T and above 1.3 K, where the upturn seen for 
T = 0.79 K in (a) is absent. 
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FIG. 14: 1/T* d versus temperature T at applied magnetic 
fields H = 1.06, 1.68, and 3.0 T for the LiV 2 4 pow- 
der samples (a) 6b with ndcfcct = 0.21 mol%, (b) 7a with 
"■defect = 0.49 mol%, and (c) 6a with n de fect = 0.83 mol%. 
The solid curves are fits to the data by Eqs. (|20[) and (1221 1, 



we assume that n = t is the same for all defects, and 
concentrate on the geometric inhomogeneity as caused 
by a varying distance between nuclear spins and defects. 
In order to evaluate the average over defect positions we 
write J^i h(Ri) = 2j <7-t/i( x i) where h(yi) is an arbitrary 
function of position y^, Xj is a vanadium site position, 
we assume the magnetic defects are at vanadium lattice 
sites, and the sum over i now runs over all vanadium 
sites, not only magnetic defect sites. The random vari- 
able gi is 1 with probability ndefect and with probability 
(1 — ndefect) and assumed uncorrelated for different sites, 
i.e. defect positions are uncorrelated. Following Refs. [27| 
andua we obtain in the limit of low defect concentration 



(ndcfcct < 1) that 
P(t) 



exp 



"defect 



1 — exp 



The sum over the lattice is evaluated as an integral. In 
the long time limit we obtain the result p (t) as given in 
Eq. © where 1 /T* d is then given bj*2L2£ 



1 _ 16tt 3 

7fT — —f— \PN "defect) 
1 ld 9 



C T 



(14) 



1 + Kr) z 

where pn is the number density of L1V2O4 formula units. 
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In the Appendix, we show that instead of solving for 
the relaxation curve, we can understand the occurrence 
of a root exponential relaxation as arising from our cal- 
culated probability distribution of nuclear 1/Ti values. 
We will discuss the temperature and field dependences 
of 1/T* d when we study the dynamics of the magnetic 
defects in Sec. IIVB5I 



3. Hole Burning Experiment and the Dependence of 
Relaxation on the Position in the Spectrum 

Bloembergen and coworkers 2 ^ have considered the 
problem of spin diffusion in the frequency domain (spec- 
tral diffusion) in a spectrum with the same kind of inho- 
mogeneous broadening as discussed for the longitudinal 
spin relaxation. The time for a hole to diffuse through 
the whole spectrum by two-spin mutual spin flip is esti- 
mated to be T^/Tlj* 3 , where T 2 is the intrinsic nuclear 
spin-spin relaxation time and T 2 * is the half width at 
half maximum of the transient echo signal. In the pow- 
der sample 6a (^defect = 0.83 mol%), T 2 « 200 p,s and 
T 2 * 5 fj,s, so T 2 4 /T 2 * 3 = 32 s. In the powder sample 6b 
(n dofoct = 0.21 mol%), T 2 « 200 ^s and T* « 20 ^s which 
give Tjf/TJ 3 = 200 ms. Both diffusion times are much 
longer than the values of T* at 4.2 K in each sample in 
Figs. [7] and [8] and are thus consistent with the lack of 
spectral diffusion in Fig. [TO] 

The higher relaxation rates at the wings of the spec- 
trum compared to that at the peak of the spectrum 
as shown in Fig. [11] can also be qualitatively explained 
within the approach where we only include the geomet- 
ric distribution of the nuclear spin to defect separations. 
For concreteness of discussion, we assume that the local 
field is purely dipolar. Denote the angle between the ap- 
plied magnetic field and the direction from a magnetic 
defect to a nuclear spin by and the distance between 
the defect and the nuclear spin by r. The NMR frequency 
shift depends on and r through (1 — 3 cos 2 6>)/r 3 , while 
the nuclear spin-lattice relaxation rate depends on 9 and 
r through sin 2 6* cos 2 9/r 6 ^- The higher relaxation rates 
observed at the wings compared to that at the peak of 
the spectrum is due to the monotonic decrease of both 
the frequency shift and the nuclear spin-lattice relaxation 
rates with increasing distance r. The nuclear spins with 
larger frequency shift will also have a higher probability 
of having larger 1/Ti values. 

4- 7 Li Nuclear Spin Diffusion 

The p(t) of powder sample 6b (refect = 0.21 mol%) 
in Fig. [T2T a) deviates from a root exponential decay at 
t > 100 ms at T = 0.61 K. This can be attributed to the 
effect of spin diffusion.— Spin diffusion tries to establish 
a common spin temperature (i.e., the same longitudinal 
magnetization) among nuclear spins at different distances 
from the defects and results in a single exponential relax- 




t(ms) 

FIG. 15: Semilog plot of the nuclear spin relaxation function 
p(t) in Eq. (|8]) versus time t after saturation for LiV204 pow- 
der sample 6b with nd e f ec t = 0.21 mol% at H — 1.06 T and 
T = 0.61 K. The straight line is a single exponential fit to 
the data at t > 110 ms. 



ation at long t. Figure [15] displays p{t) versus t of the 
same data as in Fig. H2T a) at T = 0.61 K, but on a 
semilog scale, which suggests a single exponential decay 
at t > 100 ms. A fit by p{t) = Aexp(-t/T t ) to the data 
at t > 110 ms gives 1/Ti = 1-1 s" 1 and A = 0.86. The 
best fit is shown as the straight line in Fig. [TO] 

A crossover from a root exponential to a single expo- 
nential decay occurs in the case of diffusion limited re- 
laxation as discussed first by Blumberg in Ref. The 
time t c , at which the crossover from a root exponential 
to a single exponential decay takes place, is related to 
the spin diffusion constant D through 30 

t c = C l l 2 D- 3 ' 2 , (15) 

where we dropped the subscript of Cj in Eq. (TTTj) since we 
assume the same spin dynamics for all the magnetic de- 
fects. The diffusion constant D is related to the rate W of 
mutual flips of nearest neighbor nuclear spins through 20 

D = Wa 2 , (16) 

where a is the distance between the two spins. The rate 
of the single exponential decay at long times in Fig. [15] 
is given by 3 ^ 

i = S.bpNn^C^D 3 ^, (17) 
J-i 

where pn is the number density of L1V2O4 formula units. 

In order to confirm the spin diffusion interpretation, 
below we will show that the estimated crossover time 
t c and 1/Ti are of the same order of magnitude as the 
observed t c ~ 100 ms and 1/Ti = 1.1 s^ 1 , respectively. 
The mutual spin-flip is due to nuclear dipolar interactions 
and the value of W can be estimated using Fermi's golden 
rule. For nuclear spins having / = 1/2, after averaging 
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over the angular dependence, one obtain, 



2 -y 4 h 2 



,20 



(18) 



where p(0) is the spectral density of the two spin system 
at zero Zeeman energy and j n is the gyromagnetic ratio of 
the nuclear spins. The 7 Li nuclei have spin / = 3/2, but 
an expression for W when / = 3/2 is not available, and 
the above equation for W should provide at l east a rou gh 
estimate of W. Approximating p(0) by 1/ y/2ir(Aui 2 ), 20 
where (Aw 2 ) = 288 kHz 2 is the Van Vleck second mo- 
ment of the 7 Li nucleip^ and taking a = 3.57 A, which is 
the nearest-neighbor 7 Li- 7 Li distance in L1V2O4, we find 
W = 46 s- 1 from Eq. JTSJ and D = 5.9 x lCr 14 cm 2 /s 
from Eq. (fl75|) . 

The value of C in Eq. (fT5|) can be obtained from 



Eq. (TTJ} where 1/T* d is measured using Eq. (0 from 
the initial root exponential part of p(t) in Fig. [T2f a) for 
powder sample 6b. At T = 0.61 K and H = 1.06 T, 
one obtains 1/T* d = 0.7 s" 1 , so one has C = 1.4 x 
10~ 41 cm 6 s _1 . Using the above D — 5.9 x 10~ 14 cm 2 /s 
and ndefect = 0.21 mol%, Eq. (fTS"]) yields the crossover 
time t c = 220 ms and Eq. (fT7|) yields the long time decay 
rate 1/Ti = 1.9 s _1 . Due to the uncertainty in our esti- 
mate of the parameter D and the approximate nature of 
Eq. (fTS"]) . the estimated t c and 1/Ti values should be con- 
sidered to be consistent with the observed t c ~ 100 ms 
and 1/Xi = 1.1 s -1 , respectively. 

The absence of a deviation from root exponential be- 
havior in samples 7a (n dc f cct = 0.49 mol%) and 6a 
(^defect = 0.83 mol%) as shown in Figs. [T2Tb) and (c) 
may be due to the effect of inhomogeneous broadening, 
which decreases the probability of overlap of Zeeman level 
splittings of neighboring 7 Li nuclei and results in a de- 
crease in the spin diffusion constant D. Furthermore, 
due to the higher concentrations of the defects, values of 
p(t) at t > 100 ms in these two samples are much smaller 
than in the 0.21 mol% sample, making such a deviation 
more difficult to observe. 



5. Magnetic Defect Spin Dynamics in the Powder Samples 

In this section, we discuss the relation of the nuclear 
spin-lattice relaxation rate to the dynamics of the mag- 
netic defects in the powder samples. We consider the 
weak collision limit h <C H , where h is the magnitude of 
the local fluctuating field at the nuclear site and H is the 
magnitude of the static applied field. The nuclear spin- 
lattice relaxation rate 1 jT\ due to an electronic magnetic 
defect spin at the origin is then given by 2 ^ 

1 1 f°° 

Y {r) ^¥ £ A ^ r H (S a (0)S a (t))eM^nt)dt, 

^ a—x.y,z ~°° 

(19) 

where r is the position of the nuclear spin with respect 
to the magnetic defect, A a (r) is the hyperfine coupling 



constant between the nuclear spin and the magnetic de- 
fect, uj n = 7Li-ff is the nuclear Larmor angular frequency, 
and (S a (0)S a (t)} (a = x,y,z) are the magnetic defect 
spin autocorrelation functions. 

As indicated in the ac susceptibility measurements, the 
peaks in 1/T* d versus T in Fig. [14] are related to spin 
freezing of the magnetic defects. As a first attempt, we 
assume a single exponential decay for the magnetic de- 
fect spin autocorrelation functions and assume that the 
freezing process is due to an energy barrier so that the 
correlation time r follows 



A 

T exp [ — 



(20) 



where To is the fluctuation rate of the paramagnetic de- 
fects at high temperature and A is the energy barrier in 
temperature units. For simplicity, we will assume that 
all the magnetic defect spins in a sample have the same 
correlation time r. For dipolar or RKKY interactions, 
the nuclear spin-lattice relaxation rate of a Li nucleus 
due to a nearby defect at distance r is 2 ^ 



1 



(r) 



2i? M 2 7 2 i 5(5 + l) 



5r 6 



1 



(21) 



where the angular dependence is ignored and the prefac- 
tor is written in such a way that R = 1 would correspond 
to relaxation due solely to the fluctuating dipolar field of 
the longitudinal component of the magnetic defect spin. 
The presence of additional relaxation channels would in- 
crease the value of R. By comparing Eqs. (fTO]) . (fTTj) . and 
(l2"Tj) . one has 

Cb = 1^-^5(5+1). 

Inserting this expression for Co into Eq. |[T4|). the mea- 
sured characteristic relaxation rate 1 /T* d can be written 
as 



1 



T, 



= R- 



ld 



,32tt 
~45" 



Mb7Lx 5 ( 5 + l )PN n defect 



1 



(22) 



At high temperatures, r is generally much shorter than 
the inverse of the nuclear Larmor frequency l/w n . As r 
increases with decreasing temperature T, a peak appears 
in 1~/Tf d versus T at the temperature where r = l/o; n . 

We fit the 1 /T* d data in Fig. [14] on all three powder 
samples simultaneously by the combination of Eqs. (j2"D|) 
and (f22|) . Possible field and temperature dependences of 
the parameter R are ignored in the fit. There are seven 
free parameters in the fit, R and A for each sample and tq 
which is assumed to be sample independent. The fitting 
results are displayed in Fig. [TJ] by the solid curves. The 
best fit value of tq is 4.1 x 10~ 10 s and the best fit values 
of R and A for each sample are listed in Table[TT] The en- 
ergy barrier A increases with increasing concentration of 
magnetic defects, which indicates that the dynamic slow- 
ing down with decreasing temperature originates from 
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TABLE II: Best fit values of prefactor R and energy barrier 
A obtained by fitting the 1/T* d data in Fig. [14] by a combi- 
nation of Eqs. (|20|l and (|22[) . In order to see the correlation 
between defect concentrations and R and A, the values of 
"■defect 1/5(5* + 1) are also listed. 



Sample 


n d efect\/'S'(5' + 1 ) (mol%) 


R 


A (K) 


6b 


0.85 


0.04(4) 


1.1(1) 


7a 


1.9 


0.17(6) 


1.8(2) 


6a 


3.6 


0.24(4) 


2.5(2) 



the interaction between the magnetic defects. Interaction 
between magnetic defects should increase with increasing 
concentration of the magnetic defects since the average 
nearest-neighbor distance decreases. 

The values of R in all three samples are much less 
than unity, a fact which cannot be explained by the pres- 
ence of other nuclear spin-lattice relaxation mechanisms 
since additional relaxation mechanisms would increase 
R. Such small values of R might be related to the spin- 
glass like freezing as observed in the ac magnetic suscep- 
tibility measurements. In spin glass systems, the spin 
autocorrelation functions are highly nonexponential, 31 ' 32 
which reduces the spectral density of the magnetic defect 
spin fluctuations at tj n as compared to the Lorentzian in 
Eq. (f2"2"]l . The reduction in spectral density thus results 
in a reduction in the fitted value of R in Eq. 



C. Effects of a Size Distribution of Magnetic 
Droplets 

Our previous discussion demonstrated that the non- 
exponential time dependence of the 7 Li nuclear spin- 
lattice relaxation in the sample of single crystals below 
1.3 K cannot be understood in terms of the interaction 
of the nuclear spins with separable contributions from a 
heavy Fermi liquid and a random distribution of point- 
like magnetic defects. This problem is exemplified by the 
unphysical p(t) data at long times in Fig. ll3f a) where the 
data move away from equilibrium rather than towards 
equilibrium with increasing time t. An appealing ap- 
proach that refines our previous analysis is based on the 
hypothesis that the magnetic defects are not point-like, 
but are rather magnetic droplets with an average size 
significantly larger than an atomic size. This hypothesis 
is supported by the large average spins of the magnetic 
defects inferred above in Sec. IIII Al and previously- 
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to 

be S ~ 2-4. We envision that the magnetic droplets have 
a variable size and spin and that the corresponding mi- 
croscopic internal relaxation time r varies with droplet 
size. Our arguments parallel those for the unusual spin 
dynamics due to statistically rare fluctuations that be- 
comes crucial in the context of Griffiths singularities close 
to phase transitions. 33 ' 34 

In most three-dimensional systems a high tempera- 



ture Curie- Weiss behavior of the susceptibility with Weiss 
temperature 9 is indicative for magnetic ordering at a 
lower temperature T ~ 8. However, the geometric frus- 
tration for antiferromagnetic ordering within the vana- 
dium sublattice of the spinel structure is likely the reason 
why long-range antiferromagnetic order is suppressed in 
pure LiV2 04 and a heavy electron state emerges instead. 
Defects in the crystal structure can locally lift the frus- 
tration and easily cause magnetic order in a finite region 
of volume ~ £ d around the crystal defect where the lin- 
ear size £ is, due to the proximity to an ordered state, 
expected to be larger than the interatomic spacing. Here 
d is the dimensionality of the system. These finite re- 
gions of magnetic order are what we are calling magnetic 
droplets. 

Fluctuations in the local tendency towards order usu- 
ally lead to a probability density for the linear droplet 
size £ that decreases exponentially with the volume of 
the droplet. In three dimensions one obtains^ 



p(0 



— exp 



(23) 



where £0 is the mean droplet size. It is then natural 
to assume that the typical internal excitation energy of 
a droplet decreases with increasing droplet size. These 
excitations can change the magnitude and/or direction of 
the magnetic moment of the droplet. Often the excitation 
energy e ~ H/t varies with £ according to a power law 



1 = DC 



(24) 



where we have set h = 1 and r is the relaxation time of 
the excitation. The arguments of Ref. [H yield a result 
ip = d for classical Heisenberg spins, where d — 3 for the 
present problem. Quantum effects can then yield devia- 
tions from this behavior and typically yield exponents ip 
that are larger than the above classical result^ In some 
cases, such as Ising spins in a magnetic field or Heisen- 
berg spins in a metallic hostf2£ the quantum dynamics 
can lead to a droplet dynamics in three dimensions where 



(25) 



Finally, in case of magnetic defects with Ising anisotropy 
inside a metallic host, the quantum dynamics of the de- 
fect leads to a freezing of all droplets beyond a certain 
size, typically of order £ (see Refs. H3 and|H). 

To illustrate the effects of such droplet dynamics we 
concentrate first on the case Eq. (|2"4"|) . In the absence 
of a microscopic model of the droplet spin dynamics we 
leave ip an open parameter of the model. We start from 
Eqs. © and HU) and write 



p(t) = J d£P(£)exp 



t 



Tt (o 



1/2' 



(26) 



where T-j* (£) 1 is given in Eq. (f2"2")l with r replaced by 
t (£). In the evaluation of this average we have two very 
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distinct limits. If w n r <C 1 the long time behavior of the 
system is dominated by small clusters since T^*" 1 oc r and 
slow nuclear relaxation is caused by fast droplet dynam- 
ics. In this limit, the average over droplet sizes will not 
cause any changes in the stretched exponential behavior 
as compared to our previous, purely geometric, consid- 
erations for point-like magnetic defects. The inequality 
lu u t 1 is expected to be valid at higher temperatures, 
before a freezing or dramatic slowing down of the droplet 
moments sets in with decreasing temperature. This is 
fully consistent with our findings that we obtain f3 — i 
in this regime. The situation changes dramatically in the 
limit where uj n T 3> 1, relevant at lower temperature, i.e. 
most likely for T < 1.5 K. Now T*~ l oc r" 1 and slow 
nuclear relaxation is tied to slow droplet relaxation. In 
this regime we find 



1 



T* (0 



-V 



(27) 



In the long time limit, the average over the droplet sizes 
in Eq. (|26[) can be performed via saddle point integration 
and yields 



p (t) = exp 



t 



1, droplet 



with 



= 



6 + 



(28) 



(29) 



i.e., the additional inclusion of droplet size variations 
yields a stretched exponent (3 < \. The characteristic 
relaxation rate of the droplets is 

1 ^ 327T 3 ^7gjg(^+l)p>Lect (30) 



T, 



1, droplet 
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We emphasize that the static size distribution of mag- 
netic droplets only becomes apparent by NMR at low 
temperatures such that uj n T > 1, where u n is the nuclear 
Larmor frequency. We also emphasize that the analy- 
ses in this section do not include interactions between 
droplets and a possible resultant spin glass phase. Above 
T ~ 1.5 K size variations of the droplets will not affect 
the long time nuclear relaxation, even if they are present. 

In case of an exponential dependence of the droplet 
excitation energy, Eq. (|2"5"|). the long time dynamics of 
the nuclear spins at low T is even more dramatically af- 
fected and changes to a power law decay p (t) oc t~ x , with 
nonuniversal exponent A. At the same time heat capacity 
and susceptibility measurements also experience power 
law behavior^ In this regime droplet quantum fluctu- 
ations will not only dominate the long time dynamics 
of the nuclear decay, but also thermodynamics quanti- 
ties. This might be responsible for the fact that it is not 
possible any longer to clearly separate the response of 
the underlying Fermi liquid from that of the droplet at 



low temperatures in our single crystals, as indicated by 
Fig. [Ha). 

Our data for single crystals do not allow at present 
to distinguish between the different scenarios outlined in 
this section. They do however suggest that dynamics of 
the magnetic droplets plays an important role. 



SUMMARY AND CONCLUSIONS 



The modeling approach in Sec. IIVBI gives a good de- 
scription of our 7 Li NMR results from 0.5 to 4.2 K for 
our LiV204 powder samples containing magnetic defect 
concentrations up to 0.83 mol%. This approach assumes 
that (i) there is a random distribution of magnetic point 
defects, (ii) the heavy Fermi liquid in magnetically pure 
L1V2O4 survives in samples containing up to ~ 0.8 mol% 
magnetic defects, and that (iii) the influences of the mag- 
netic defects and of the Fermi liquid on the magnetization 
and NMR properties are separable. This description ex- 
plains very well the defect concentration-independent \o 
value from our low-temperature magnetization measure- 
ments, the inhomogeneous broadening of the 7 Li NMR 
spectrum, the nonexponential 7 Li nuclear spin-lattice re- 
laxation versus time behavior, and the lack of spectral 
diffusion in the 7 Li NMR hole burning experiments. It 
also explains the smaller 7 Li nuclear spin-lattice relax- 
ation rate at the peak of the spectrum as compared to 
that at the wings. However, it is hard to reconcile the 
picture of magnetic point defects with the high magnetic 
moments for the defects (spins of 2-4) deduced here (see 
Tabled and in Refs. and[T3 from magnetization mea- 
surements. These large defect spin values suggest that 
the magnetic defects may not behave like point-like mag- 
netic moments under all circumstances. In Sec. HVCl we 
discussed magnetic defects that are more extended enti- 
ties that we have called magnetic droplets with a distri- 
bution of sizes, and likely a distribution of spin values. 
We showed that such a size distribution can affect the 
NMR magnetization recovery at long times at low tem- 
peratures in the regime where interactions between the 
magnetic droplets can be neglected. 

Our study shows that there can be different kinds of 
magnetic defects in the L1V2O4 system. As revealed by 
the nuclear spin lattice relaxation rate data and ac mag- 
netic susceptibility measurements at 14 MHz, it is amaz- 
ing that the magnetic defects in the powder samples un- 
dergo a spin glass-like freezing below 1 K, while the mag- 
netic defects in the single crystals with a similar magnetic 
defect concentration exhibit a very different behavior at 
such low temperatures, with no evidence for spin freez- 
ing. The different kinds of magnetic defects and/or inter- 
actions in the crystals and powders must be associated 
with different types of structural defects in the system, 
which might be expected because the crystals are grown 
at about 1000 °C whereas the powders are synthesized 
at 700° C. Different types of magnetic defects were even 
found in an annealing study of different single crystals^ 
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where heat treatment at 700 °C was found to remove the 
magnetic defects in one but not in other single crystals. 

Further experiments on the single crystals arc urgently 
needed at ~ 1 K and below. These experiments should 
test whether the Fermi liquid is modified by quantum 
fluctuations of large magnetic droplets in the single crys- 
tals at T < 1.3 K, whether the magnetic properties of the 
crystals contain separable contributions from the Fermi 
liquid and the magnetic defects, and whether the conduc- 
tion electrons in the crystals even form a Fermi liquid. In 
addition, the origin of the distinct kinks at about 1.4 K 
in the temperature dependences of 1/T-j* and f3 in Fig. [7] 
and of M(oo)T in Fig. [5] for the single crystals remains 
to be explained. 
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APPENDIX: PROBABILITY DISTRIBUTION 
UNDERLYING STRETCHED EXPONENTIAL 
RELAXATION WITH = 1/2 

Here we discuss the stretched exponential 7 Li nuclear 
relaxation versus time following saturation that arises 
from interactions between the 7 Li nuclear spins and a 
low concentration of magnetic defect spins, which is p(t) 
in Eq. ©. We demonstrate that instead of theoretically 
solving for the relaxation curve and showing that it is a 
stretched exponential with (3 — 1/2 as in Sec. IIV B~2| we 
can understand the occurrence of this root exponential 
relaxation as arising from the probability distribution of 
nuclear 1/Ti values. 

We assume that p(t) in Eq. © is due to a contin- 
uous sum of exponential decays with a distribution of 
relaxation rates 1 jT\ . Then one can write the stretched 
exponential relaxation function in Eq. ^ as 

/>oo 

e~W T ^= P( S ,P)e- st / T ^d S , (A.l) 



where s equals 1/Ti normalized by 1/T* d , i.e., s = 
T£ d /Ti, and P(s,(3) is the probability density for occur- 
rence of s for a fixed exponent (3 with < f3 < 1. Thus 
the stretched exponential function is the Laplace trans- 
form of P(s, (3). Closed analytic expressions for P(s, (3) 
with rational values of (3 can be obtained from the inverse 
Laplace transform of the stretched exponential function, 
and physical interpretations of the parameters 1 /T* d and 
(3 have been determined * 12 ' 26 We show below that the 
probability distribution of 1/Ti due to dipolar interac- 
tion of nuclear spins with dilute magnetic defects corre- 
sponds very well to the 1/Ti distribution leading to the 



stretched exponential relaxation in Eqs. ([9]) and (|A. 1[) 
with 13 = 1/2. This probability density isi 2 -^ 



P(s,l/2) 



e "a 



Airs 3 / 2 



(A.2) 



A plot of this distribution function is given below as the 
solid curve in Fig.[T?5J P(s, 1/2) is proportional to s~ 3 / 2 
for large s [in general P(s,/3) at large s is proportional 
to s _ ( 1+ ^], and has a low-s cutoff since e _1 / 4s expo- 
nentially approaches zero at small s values (also true in 
general for arbitrary /3)i2S 

A qualitative s -3 / 2 dependence of the 1/Ti probability 
distribution arises due to a r~ 6 dependence of I/Tl as in 
Eq. (JTUJ) as follows. Ignoring the angular dependence in 
Eq. ([TO]) and assuming the same dynamics for all the 
magnetic defect spins such that Ci = C, in the single 
paramagnetic center limit the distribution of s arising 
from a continuum distribution of nuclear spins around a 
magnetic defect ia~2. 



Pgoo(s) cx r 2 ^ 
as 



CT* \ 1/6 oc s 



-3/2 



(A.3) 



This s dependence is the same as the large-s limit of 
the result in Eq. (|A.2[) for stretching exponent (3=1/2 
noted above. The distribution (|A.3|) diverges as 1/Ti ap- 
proaches zero. This divergence is caused by the single 
impurity approximation. Nuclei with 1/Ti approaching 
zero correspond to those far away from the paramagnetic 
center. Due to the finite distance between different para- 
magnetic centers in a system with a finite concentration 
of them, the actual probability of finding a nuclear spin 
with 1/Ti — * should instead vanish, so a low 1/Ti cutoff 
has to be applied. This results in a distribution function 
with a shape (see the dashed curve in Fig. [TTJ] below) 
roughly similar to that in Eq. (|A.2[) (the solid curve in 
Fig. [TJ below). 

We have carried out a numerical simulation of the 1/Ti 
probability distribution that turns out to be in very good 
agreement with the exact result for the probability distri- 
bution in Eq. (|A.2[) for the stretched exponential function 
with f3 = 1/2. In the simulation, we calculated the 1/Ti 
distribution of 7 Li nuclei due to a random distribution 
of dilute point-like paramagnetic centers (defects) in the 
L1V2O4 spinel structure. The paramagnetic defects ran- 
domly occupy the vanadium sites with a probability of 
0.25% (rtdefoct = 0.5 mol%) and the configuration of the 
random defects repeats every 80 unit cells in all three 
crystallographic axis directions. The 1/Ti of each 7 Li 
nucleus is calculated using 2 ^ 



1 

Ti 



C 



15 sin 2 9i cos 2 9, 1 



(A.4) 



where C is a constant, rj is the distance between para- 
magnetic center i and the 7 Li nucleus, and Qi is the an- 
gle between the applied magnetic field and the vector 
from paramagnetic center i to the Li nucleus. The ap- 
plied magnetic field was arbitrarily chosen to be along 
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FIG. 16: 7 Li nuclear spin- lattice relaxation rate X/Ty proba- 
bility distributions P(l/Ti), normalized by T\, due to dilute 
point-like paramagnetic defects. The filled circles are results 
obtained from computer simulation of LiV2C>4 where the mag- 
netic defects randomly occupy vanadium sites with probabil- 
ity 0.25%. The solid curve is the best fit of these data by 
Eq. (g]H with 1/Ti* = 0.067(1) s" 1 . The dotted curve is a 
plot of p(x) = 0.8x~ 3 ' 2 with a small x cutoff at x c = 0.0256, 
where x = oq/CTi. The lower cutoff x c is chosen so that the 
normalization condition J°° P(x) dx — 1 is satisfied. 

the (001) direction. Equation (|A.4[1 is the nuclear spin- 
lattice relaxation due to the dipolar magnetic field fluc- 
tuations from the longitudinal spin component of the 
paramagnetic defects^ In the presence of a strong ap- 



plied magnetic field, the transverse spin fluctuation is 
often modulated by the Larmor frequency of the elec- 
tronic spins and thus has negligible contribution to 1/Ti 
for which magnetic field fluctuations at the nuclear Lar- 
mor frequency are most important^ The summation 
over i in Eq. (|A.4p includes all magnetic defects with 
ri < 20ao, where ao = 8.24 A is the lattice constant of 
cubic L1V2O4. We checked that changing the summation 
range to r,; < 10ao resulted in a negligible change in the 
1/Ti distribution. 

The distribution of 1 jT\ resulting from the above sim- 
ulation is displayed as the filled circles in Fig. [TBI where 
the distribution is normalized by T\. The maximum re- 
laxation rate plotted in Fig. [TB] at a^/CTi w 20 is not 
a large relaxation rate cutoff to the probability distri- 
bution. Data at larger relaxation rates are not plotted 
due to insufficient statistics. The simulated 1/Ti distri- 
bution in Fig. [TH can be fitted very well by Eq. (|A.2[) 
with 1/Ti* d = 0.067(l)C/ag, as shown by the solid curve 
in Fig. [TH1 1/T* d calculated from Eq. (|14p is equal to 
0.088C/ag, close to the simulated result. The difference 
may be due to the neglected angular dependence in de- 
riving Eq. (|14p . For comparison, a simple power law dis- 
tribution P(x = al/CTi) = 0.08a;- 3/2 with a small x 
cutoff of x c = 0.0256 is also displayed in Fig. [TCJ This x 
dependence is the same as in Eq. (| A.3|) and is the same 
as the asymptotic large-s dependence of Eq. (|A.2[) . The 
prefactor 0.08 is chosen to make the distribution over- 
lap with the simulated result at large x, and the low-x 
cutoff x c — 0.0256 is determined from the normalization 
condition J^ 00 P(x) dx = 1. 
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